home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 1999 August / SGI Freeware 1999 August.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / bead.ode.z / bead.ode
Encoding:
Text File  |  1998-10-28  |  535 b   |  24 lines
  1. # This example simulates a bead sliding on a smooth circular wire.  The
  2. # output displays the height of the bead as a function of time.
  3. #
  4. # You may run this example by doing:
  5. #
  6. #    ode < bead.ode | graph -T X -C
  7. # 
  8. # or alternatively, to get a real-time plot,
  9. #
  10. #    ode < bead.ode | graph -T X -C -x 0 5 -y 0 2.5
  11. #
  12.  
  13. a = 1  # radius of circular wire
  14. g = 10 # acceleration due to gravity
  15. w = 10 # angular velocity of circular wire
  16.  
  17. the' = vthe
  18. vthe' = (w^2)*sin(the)*cos(the) - (g/a)*sin(the)
  19. the = 0.1
  20. vthe = 0
  21.  
  22. print t, the
  23. step 0,5
  24.